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We consider anisotropies in tiie stochastic background of gravitational-waves (SBGW) arising 
from random fluctuations in the number of gravitational-wave sources. We first develop the general 
formalism which can be applied to different cosmological or astrophysical scenarios. We then apply 
this formalism to calculate the anisotropies of SBGW associated with the fluctuations in the number 
of cosmic string loops, considering both cosmic string cusps and kinks. We calculate the anisotropies 
as a function of angle and frequency. 
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A stochastic background of gravitational-wave 
(SBGW) radiation is produced by a large number of 
weak, independent and unresolved gravitational wave 
sources. The sources of the SBGW can be isotropic 
or anisotropic. For the case of sources of cosmological 
origin [l|-|3[ the distribution of the gravitational wave 
sources is expected to be isotropic, while astrophysical 
sources such as rotating neutron stars 0] or magnetars 
may have an anisotropic distribution. Even in the 
case of an a priori isotropic source distribution, random 
fluctuations in the number of sources will (in general) 
give rise to anisotropies. Such anisotropies are analogous 
to the anisotropies observed in the cosmic microwave 
background (CMB) radiation and would carry additional 
information about the gravitational-wave sources that 
generated them. 

Cosmic strings are expected to contribute to the SBGW. 
Cosmic strings are predicted by a large class of unified 
theories ^Qr^ as remnants of spontaneously broken 
symmetries at phase transitions in the early Universe, 
as well as in string-theory-inspired cosmological scenar- 
ios [§] . Once formed, a network of cosmic strings evolves 
toward an attractor solution called the scaling regime, 
in which the statistical properties of the network, such 
as the average distance between strings and the size of 
loops at formation, scale with the cosmic time. The 
gravitational interaction of strings is characterized by 
the dimensionless parameter Gfj,, where G is Newton's 
constant and /z is the tension. The current CMB bound 



on the tension is G/i < 6.1 x 10~^ [l3, III- Cosmic 
string cusps, regions of string that acquire enormous 
Lorentz boosts, are expected to generate large transient 
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gravitational-wave signals 
bursts could be observable by current and planned 
gravitational wave detectors [15( for values of G/i as low 
as 10"^'^, which may provide a probe of a certain class of 
string theories @ . A SBGW produced by the incoherent 
superposition of cusp bursts from a network of cosmic 
strings and superstrings was considered in [l^, [l7l| , and 
it was later shown that kinks contribute to the SBGW 



at the same order as cusps [18|. These sources of SBGW 
are also observable by current and planned detectors, 
for a wide range of the parameter space, see 17 , 3] and 
references therein. 

In this paper, we develop a general formalism to treat 
SBGW anisotropies. In particular, we consider two-point 
correlations in SBGW between two different directions 
in the sky, which arise from random fluctuations in 
the number of gravitational-wave sources. While this 
formalism is applicable to a variety of cosmological and 
astrophysical SBGW models (see [19[ and the references 
therein), we illustrate it for the specific case of cosmic 
(super)string cusps and kinks. 

Anisotropies in the SBGW: We start from formalism 
in references [13, 2A\ and extend it to treat angular 
dependence. The energy density of SBGW at frequency 
/ corresponding to sources in the direction Cl is given by 
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where dpgw is the energy density of gravitational waves 
in the frequency range f to f + df and pc is the criti- 
cal energy density of the Universe. Let us assume that 
sources are characterized by a set of parameters ^ - in the 
case of cosmic strings, redshift z is one such parameter. 
Therefore ilgw is an integral over the parameter space C 
which we propose to discretize as follows: 

C!gw(/,f2) = J dCnic,n)wifX,^) 

i 

= J2NiC^^^MfX^:^)■ (2) 

i 

Here we assume that the parameter space can be divided 
into disjoint volumes A(^i), centered at Ci: whose size 
is large compared to the correlation length of the num- 
ber of sources. In other words, a statistical fluctuation 
in the number of sources in one volume would have no 



implications on the number of sources in any other vol- 
ume. We further define n(Ci, ^) as the number density of 
sources (i.e number per parameter space volume) in the 
direction Q with the parameter set Ci, and w{f,Ci,^) as 
the contribution to Jlgw of one source at frequency /, 
in the direction Cl, and with the parameter set Ci- We 
also define A^(Ci, ^) = n{<^i,^l)A{<^i) as the total number 
of sources with the parameters in the range from to 
Ci + A{Ci) and in the direction fl. The contribution of 
one source is given by 
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where h{f,Ci,D,) is the strain of the gravitational wave 
with frequency / originating from a source with param- 
eters Ci and at the line of sight Ct. TZ{f, Ci, ^) represents 
the observable part of the gravitational radiation from 
the source, i.e. it incorporates the propagation of the 
wave in the expanding universe as well as possible beam- 
ing effects, see [l2. 13|. 

The angular dependence of N can originate from 
anisotropic source distribution. Moreover, even in the 
case of an a priori isotropic SBGW, random fluctuations 
in the number of sources will (in general) give rise to 
anisotropics. Note that N{Ci,^) are dimensionless num- 
bers, which are by construction uncorrelated for differ- 
ent values of the index i. Assuming Poisson distribu- 
tion, the statistical fluctuations of N{Ci,(l) are of order 

N{Ci,^)- The corresponding fluctuation in ftg^ is 



The two-point correlation of Sflg^if, i^) at two different 
directions reads 



c = (^6ng^if,n)5ng^if,n')j (5) 

= E ^(/' C^,^Mf, C,,tt') {SN{C^,(l)SN{C,,^') 



Since the fluctuations in the number of gravitational- 
wave sources are Poissonian, we propose the following 
bilinear expectation: 



6N{C^, n)6N{Q,n') ) ^ N{C^, fi) J-(7, C^) h 
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where 7 is the angle between VL and Cl' . is a function 
that incorporates the correlation properties of the grav- 
itational wave sources. Although the precise form of F 
will depend on the problem at hand, we can discuss sev- 
eral properties of this function. Firstly, we expect to see 
the maximum correlation if the two sources are close to 
each other in the physical space as well as the parameter 
space. Therefore J- must assume its maximum value at 



7 = 0, and it should decrease for larger values of 7. Since 
J- constrains 7, the angle between $7 and il', to small 
values, we keep only Cl at the right hand side of Eq. 
This is a good approximation as long as N changes slowly 
with Cl. Below we will consider an explicit example and 
discuss form of J- in more detail. Inserting this into Eq. 
(0) gives 



c = J2MCMC.,^W{f,C.,^)J'{i,C. 



dCn{c,ci)w^{f,cMmiX), 
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where we take the integral limit of the sum. This is 
a general expression applicable to both cosmological 
and astrophysical problems in which the correlation 
properties of the sources are specified by the function J". 

Cosmic Strings Case: We now apply this formalism to 
the case of cosmic strings, including gravitational- wave 
bursts from cusps and kinks, in which the distribution of 
sources is specified by the redshift z. It is convenient to 
parameterize physical quantities in terms of redshift. To 
this end, we define the following dimensionless cosmolog- 
ical functions: 
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where n{z) = A/ri„(l -I- z)^ + flril + z)* + is the 
Hubble function with Qm = 0.25, il/? = 4.6 x 10^^ and 
r^A = 1 — — ^R- We now explicitly construct the 
integral in Eq. ([7]). Firstly, the parameter space volume 
dC in this case is simply the co-moving volume. It can be 
written as HQ^ipv{z)dz, where Hq is the present value 
of the Hubble constant. This converts the co-moving 
differential volume r^dr to the corresponding differential 
volume as a function of the redshift. The next quantity 
in Eq. ^ is the number density of the loops. If the loop 
size is determined by gravitation back- reaction 12, llSj , 
the loop number density is given by 
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c{z) 
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where p is the reconnection probability, F = 50 is a di- 
mensionless parameter proportional to the power emitted 
in gravitational waves by cosmic string loops, and t{z) is 
the cosmic time which can be written as t{z) = ipt{z)/Ho. 
The function c(z) = 1 -|- (zgg — 5440) accounts for 

the fact that the loop density in radiation domination is 
about 10 times that of the matter domination. In or- 
der to define the J^- function in Eq. ([7]), we assume that 
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FIG. 1: Top: Normalized Correlation A/'C as functions of j/n for cusps and kinks for / = IHz, / = lOHz and / = lOOHz, for 
G/i = 1.0 X 10"^ p = 1 and e = 1.0 X 10"". 

Bottom: A/'C as functions of frequency for cusps and kinks, for various values of 7 and for the same model parameters as above. 
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FIG. 2: MC for cusps and kinks in Log{e) (vertical axis) Lo(/(G/i) (horizontal axis) parameter space at frequencies applicable 
to ground-based detectors (10 Hz) |l^, satellite-based detectors (1 mHz) [23|, and pulsar-based observations (10~* Hz) [2^ . 
The (base 10 logarithm of) numerical values of J\fC are denoted in the color bar for each plot. 



the number density of cosmic string cusps and kinks at 
a given z is correlated over the length scale R{z) given 
by the Hubble size, R{z) w t{z) [23|. The angular size 
spanned by this length scale at the distance r{z) can be 
calculated using the standard angular diameter-redshift 
relation as 



7z = 2 arctan 



(l + z)i?(z) 



r{z) 



(10) 



where r(z) can be written as r{z) — (pr{z)/ Hq. There- 
fore, for the given redshift z, two directions on the sky 
are correlated if their angular separation, 7, is less than 



72 . This condition can be imposed by the J- function, 
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which vanishes if 7 is larger than 7^ , the angle subtended 
by the length scale R{z). We emphasize that the correla- 
tions considered here are large scale, and arise from the 
fluctuations in the number of cosmic string loops in an 
evolving cosmic string network. This is different from 
the correlations associated with the correlation length of 
a single cosmic string loop, which are important in deter- 
mining the cosmic string signatures in the CMB [20[ . For 
cusps and kinks on cosmic string loops with sizes given 
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by the gravitational back-reaction scale we have [IE 



Wc{f,z) = 
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where a = eTG/j, is the parameter that sets the length of 
the loops. Since the function w has no angle dependence 
for kinks and cups, Eq. ^ simplifies to 

C = C(/,7)- / dzH^^ipv{z)n{z)w\f,z)Tij,z) 



dzipv {z) '3, if,z)e 
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which is a function of the opening angle, 7, and the fre- 
quency only. It is important to note that large rare events 
which occur at rates smaller than the relevant time-scale 
of the experiment are excluded from C in numerical 
calculation . This exclusion removes the a priori diver- 
gence of the integrand at z = 0. The integrand of Eq. 

quickly vanishes with increasing redshift, implying 
that the dominant contribution comes from low redshifts. 
The small values of redshift correspond to closer sources, 
which have larger angular size in the sky. Therefore the 
angular dependence of correlations will be rather flat for 
small angles, and it will rapidly vanish for large angles, 
for which small values of redshift are excluded from the 
integral by 7^(7, z). In order to understand the relative 
strength of the fluctuations at a given frequency / com- 
pared to rigw(/) (integrated over all sky) we define the 

following quantity: A/'C(/, 7) = which we refer 

to as the normalized correlation. We numerically eval- 
uate the integrals in Eq. (jl3p for kinks and cusps and 
calculate the normalized correlations, as depicted in Fig. 
1. We also do a parameter scan in e — Gji space. Fig. 2 
shows the density plot for the strength of the background 
VLgyfj and MC at various values of / for cusps and kinks. 
Conclusions: In this paper we have developed the for- 
malism for calculating the spatial anisotropics in the 
stochastic background of gravitational waves associated 
with the random fluctuations in the number of sources. 
The formalism is applicable to a variety of cosmological 
and astrophysical models. We applied it to the case of 
SBGW due to cosmic (super)string cusps and kinks, and 
observed that the relative strength of the anisotropics, 
gw , can be estimated by \/\/N — \/TG]I, which 
can be as high as lO^'^. While observation of these spatial 
anisotropics is unlikely for the second-generation detec- 
tors that are currently being built (Advanced LIGO and 
Advanced Virgo), the planned third-generation detectors 



(such as Einstein Telescope) should be sufBciently sensi- 
tive to measure them over a large part of the parameter 
space. We emphasize that the general formalism devel- 
oped here can be used to distinguish between different 
SBGW models - that is, between models that predict sim- 
ilar frequency spectra and different spatial anisotropics. 
This technique will be crucial for the identification of the 
source of SBGW which is expected to be observed by the 
future generations of the gravitational- wave detectors. 
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